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Abstract 

A discussion of the stability of self gravitating Skyrmions, with a large winding number N, 
' in a Schwarzschild type of metric, is presented for the case where an isospin chemical potential 

in ' 

CN ■ is introduced. It turns out that the chemical potential stabilizes the behavior of the Skyrmion 

discussed previously in the literature. This analysis is carried on in the framework of a variational 



approach using different ansaetze for the radial profile of the Skyrmion. We found a divergent 



^ behavior for the size of the Skyrmion, associated to a certain critical value u c of the chemical 

&: 

potential. At this point, the mass of the Skyrmion vanishes. fi c is essentialy independent of 
gravitating effects. The stability of a large N skyrmion against decays into single particles is also 
discussed. 
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The stability of a Skyrme topological soliton coupled to gravity has been analyzed many 
times in the literature 111 DO . The purpose of this letter is to extend this discussion to the 
scenario where Isospin chemical potential is taken into account. 

Gravitating Skyrmions turn out to have a very rich structure, with a nontrivial spectrum 
for the asymptotic states j^, and provide a counter example to the "non hair theorem" in 
the construction of black hole solutions [3| . It is therefore interesting and natural to consider 
such objects now in the presence of chemical potentials associated to conserved charges. 

Recently, we have discussed Skyrmions as a model for baryons in the presence of isospin 
chemical potential (//), showing the existence of a critical value /i c where the mass of the 
Skyrmion vanishes, signalizing the occurrence of phase transitions. We found that the radial 
profile of the Skyrmion becomes broader asu approaches the critical value. The nucleon 
spectra was also considered in this context m. 

n 

Here we present, following the article by Glendenning et al a variational analysis for 
the stability of a self gravitating SU(2) Skyrmion for large values of the winding number 
N, showing that fi plays an stabilizing role. In particular, the parameter associated to the 
radius of the Skyrmion profile starts to grow, whereas the mass diminishes as function of 
jj,. For a certain value /i c , the radius diverges and the mass vanishes in the same way as it 
happens in the non gravitating case jj| . We found strong numerical evidence that the value 
of n c is independent of the parameters associated to gravity. 

It was found in [jj, that the Skyrmion mass behaves like 0(N 2 ). Therefore gravity 
becomes important in the large N limit. Actually, the scale where gravity becomes relevant 
is given by Nf n ~ Mpi anc k. The Skyrmion mass depends on the quotient Rs/R ~ G(Nf n ) 2 , 
where R is a measure of the Skyrmion size in the absence of gravity, whereas Rs is essentially 
the Schwarzschild radius, i.e. the natural length scale associated to gravity effects. In the 
variational approach, there is no solution for values of Rs/R bigger than a certain limit. If 
fi 7^ 0, this limit becomes bigger, signalizing the stabilizing role of the chemical potential. 
If u — > u c , we have no bound for N. 

n 

The authors of [1J found also that their solution actually does not satisfy the condition 
M(N) < NM(1), indicating that the solution with high winding number should decay into 
many Skyrmion states with lower N. In our case, the chemical potential also provides a 
source of stability near fi c . 

We would like to remark that in our previous analysis the variational approach was 
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not appropriate due to the lack of stability for large values of \i. Therefore, a full numer- 
ical treatment was unavoidable. However, in the case of the self gravitating Skyrmion, a 
variational approach is possible for any value of \i < fi c . The idea is to consider a certain 
given radial profile, which depends on a free parameter, minimizing then the mass of the 
Skyrmion. 

Let start our analysis with the nonlinear sigma-type Skyrme lagrangian 



f 2 

+ (1) 



where 



D, = d, - i| [a 3 , U] Spja, (2) 



is the covariant derivative that introduces the chemical potential p, |6| , and 



U = (T cos F(r) + id ■ n sin Fir), (3) 

is the standard Hedgehog ansatz, a are the Pauli matrices, <To is the identity and Fir) is 
the radial profile. The gravitating version (Einstein-Skyrme lagrangian) can be obtained 
through the action 

S = J (L g + L rn ) d 4 x, (4) 

with L m given by (JTJ) and 

where R is the Ricci scalar curvature and G the Newton constant. 
The explicit form of L2 and L4, without is 

= Ifl^F'irf + g^^Fir) 

+/ 3 sin 2 F(r)sin 2 #), (6) 



U,o = -^(^V 2 ^(r) 2 sin 2 F(r) 

11„33 T?t(„\2 ■ 2 77/„\ • 2 



+ ^y^F>) 2 sm 2 F(r)sin 2 
+/V 3 sin 4 F(r)sin 2 ^. (7) 

When we turn on the chemical potential, /i-dependent terms appear, so that, in an obvious 
notation, 



L m — L 2 fl + L 4:0 + L 2jIX + L itx . (8) 



with 



^ = ^.Vsin 2 F(r)sin 2 ^, (9) 



and 



= -^ 2 sin 2 F(r)sin 2 ^(^F'(r) 2 

+# 22 sin 2 F(r)). (10) 

The Schwarzschild metric depends on two arbitrary functions v{r) and A(r), in such a 
way that 



I, »» = L. 

r 2 ' r 2 sin 2 6 1 



Our starting point is the Einstein equation 

= SttGT^. (12) 
The matter energy-momentum tensor is given by 

T^ = -ri m + 2p. (13) 
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For our variational approach, where different ansaetze for the radial profile F(r) will be 
used, only the T 00 component is needed 



, 00 _ U 2 ( 2 sin 2 F(r) , F'\r) 



2 V r 2 e 2A M 
1 /sin 4 F(r) sin 2 F(r)F' 2 (r) 



2e 2 \ r 4 e 2A(r) r 2 

f 2 u 2 

JT ■ 2 771/ \ • 2 /) 

sin ^ (r J sm # 



2 

_fi 2 /sin 2 F(r)sin 2 #F /2 (r) 
~2^ V e 21 ^ 

sin 4 F (r) sin 2 #\ . 

+ -r J • ( 14 ) 

For our purposes, we will only consider the angular average over the unit sphere. Since we 
are following a variational approach, the reasons that justify the self gravitating skyrmion 
stability will be not affected by taking this angular average. In this way, we replace sin 2 (0) — > 
2/3, obtaining an averaged energy- momentum tensor component (T 00 ) 

The first component of (jl2j) reads 

r 2 G° Q = e- 2A (l - 2rV) - 1 = 87rGV 2 (T °), (15) 
where the right hand term, written in terms of the dimensionless variable x = ef n r, becomes 



87rGr 2 (T °> = AirGfl 



sin Fix) + 



sin 4 F{x) 



2x 2 



2e 2 f 2 K ; 3 

Lh 9 9 , ,2 

-x sm F(x)~ 



2e 2 /, 



2v(x) 



X 



X 



x 2 sm 2 F(x)-)F' 2 



2e 2 f 2 



+ sin 2 F(x) 



Following y|, equation (JT5Jl can be written as 



(16) 



-2A(r) 



)] = -8nGr 2 (T°) 



(17) 



A formal integration of the previous equation leads us to the identification 



r(l - e- 2A «) = 8vrG / (T°)r ,2 dr'. 

Jo 

It is natural to define the mass inside the radius r as 



(18) 



M(r) = 4tt / (T °)r /2 dr' 



(19) 



Then, we can identify the left term in (|15|) as r(l — e~ 2X ^) 2GM(r). 

In order to compare with Q], we introduce R = VS/ef^ as a relevant length scale parameter. 

In terms of R, we define the dimensionless quantities 



R 3 eu 



(20) 



In what follows we will use the dimensionless variable M(x) = ef. K GM{x). In this way, 
we find 



dM(x) 
dx 



2M(x) 
x 



p(x) + q(x), 



(21) 



where we have introduced the functions p(x) and q(x) according to 



P( x ) = 7~ET 



lRs_ ( F'(x) 
A R \ nN 
2jl 2 x 2 sin 2 F(x 



(6sin 2 F(x) +3x 2 



(22) 



and 



2/iV 2fi 2 x 2 sin 2 F(x) . 



(23) 



Different ansaetze may be used for the radial profile, in order to proceed with the varia- 
tional approach, fulfilling the boundary conditions F(0) = Ntt and F(x — > oo) = 0. Never- 
theless, the results seem to be more or less independent of the details of the parametrization 
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employed. A simple and natural profile emerges from the discussion of instanton holonomies 

a 



F(x) 



Ntt 



1 - 



X 



This profile is also valid for curved spaces 



been proved in 



. (24) 

Vx 2 + X 2 \ 

. Actually the validity of this ansatz has 



9|. Other approaches, as for example, exponential decaying profiles F(x) = 
Ntc exp(— x/X) have been used. The variational procedure leads us to the optimal value of 
the free parameter X in each case. 

For our estimates, we will employ the first ansatz. The main reason for this, is a better 
numerical convergence of the integrals which give us the mass of the Skyrmion. We would 
like to remark that this is possible when gravity effects are taken into account. In our 
previous works on the stability of hadrons described as Skyrmions and their excitations, 
in the presence of /i, this parametrization was only valid for small values of /i. In the 
present case, we do not have such instability problems when going to higher values of \i. As 
we claimed in the introduction, gravity compensates the instability of the non gravitating 
Skyrmion. 

In the large N limit, due to the rapid oscilating factor, we may approximate sin 2 F(x) ~ 
1/2. On the other side, since q(x) oc 0{1/N 2 ) and p(x) oc 0(1), we can neglect q(x), getting 



V 



dM{x) 
dx 



3 R s ( F(x 

4 ~R 



1 - 



2M(x) 



x 



P- 




(25) 



(26) 



Gledenning et al [jj found stable solutions for the Skyrmion only for values of Rs/R < 
0.26. This means that for bigger values of Rs/R, the mass of the Skyrmion does not have 
a minimum with respect to X. 

In our case, for each value of Rs/R, we are able to find a stabilizing value of /i (us), 
such that the mass of the Skyrmion has a minimum as function of X. The behavior of /j,s 
as function of Rs/R is shown in figure (JTJ). Notice that us asymptotically goes to a certain 
critical value fi c ~ 0.3/ R (80 MeV), where the size of the Skyrmion diverges. This is shown 
in figure (0), where the divergent behavior of the parameter X is plotted as function of 
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/x, signalizing also a divergent behavior for the radial extension of the profile, since bigger 
values of X imply a broader radial profile. 

In the same way, simultaneously with these effects, it turns out that M = ef n M dimin- 
ishes as function of /x, vanishing at the critical value /x c . This is shown in figure (j3J). It 
is interesting to remark that the same kind of behavior for the mass of the Skyrmion was 
found in the non gravitating scenario discussed in our previous work. This means that the 
phase transition induced by the chemical potential is triggered by the strong dynamics. 

Finally, we would like to address the problem of stability of the large N Skyrmion into 
particle emission. The condition for that is M(N) < NM(1). This means that 

dM(N) A 

d{R s /R) N' 1 1 

From figure (j3J) we see that the slope of M as function of Rs/R diminishes when /x increases. 
For ft, —>■ ft c , the slope vanishes. This means that for any value of Rs/R we can always find 
a value for /x, such that the condition (J27|) is fulfilled, implying the stability against decay. 

In this letter we have discussed stability conditions for self gravitating large N SU(2) 
skyrmions in the presence of isospin chemical potential. Our main results show that this 
chemical potential plays an important role stabilizing the solution, extending the range of 
validity for the solutions as function of Rs/R- We showed that the skyrmion mass vanishes 
at the same critical chemical potential where the size of the Skyrmion diverges. Finally, 
this scenario allows the existence of stable large N Skyrmions against decays into single 
Skyrmions. 
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FIG. 1: The behavior of fig = HsR is shown as function of R s /R. It goes asymptotically to the 
critical value jl s — > fl c . 




FIG. 2: This figure shows the evolution of the variational parameter X that minimizes the mass 
as function of the chemical potential. It diverges when \x — > fl c . 
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FIG. 3: Here we show the Skyrmion mass evolution as function of jl. The mass vanishes when 
jl -> jl c . 
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FIG. 4: Here we show the Skyrmion mass dependence as function of Rs/R for three different 
values of jl: continuous line jl = 2.5; dashed line (1 = 2.7; dash-dotted line fx = 2.85. 
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